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Abstract 

A semiphenomenological approach to the nucleoli self-energy in nu- 
clear matter at finite temperatures is followed. It combines elements of 
Thermo Field Dynamics for the treatment of finite temperature with 
a model for the self-energy, which evaluates the second order diagrams 
taking the needed dynamics of the NN interaction from experiment. 
The approach proved to be accurate at zero temperature to reproduce 
JraS and other properties of nucleons in matter. In the present case 
we apply it to determine /mS at finite temperatures. An effective NN 
cross section is deduced which can be easily used in analyses of heavy 
ion reactions. 
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1 Introduction 



The imaginary part of the nucleon self-energy £ has been the subject of intense 
research in the past |l|]. However, these studies have been done at zero temper- 
ature, of relevance for physical processes involving nucleon-nucleus scattering. 
On the other hand, in heavy ion collisions, through multiple collisions of the 
nucleons, one reaches conditions roughly similar to those of a thermal bath at 
finite temperature 0. The imaginary part of the nucleon self-energy in such a 
bath is an important magnitude which governs processes of nucleon emission, 
particle production and, in principle, most of the nuclear processes occurring 
during heavy ion collisions, including the rate of thermalization. 

The many body field theoretical treatment at finite temperature becomes, 
however, technically more involved than at T = 0. At the root of it lies the 
fact that particle annihilation operators do not give a vanishing result when 
applied to the ground state of the system, which is now a thermal distribution. 
This does not allow one to use the Wick theorem to generate the Feynman 
diagrammatic perturbation expansion, as we know for many body systems at 
T = 0, or in ordinary Quantum Field Theory in the absence of a medium. 
Yet, even if more complicated, the methods to deal with it are available. The 
traditional approach has been the one of the imaginary time formalism 0, |J . 
although lately the new method in real time formalism known as Thermo Field 
Dynamics [||, ^| is proving quite efficient and becomes more and more widely 
used. A recent review on Thermo Field Dynamics with some applications to 
nuclear matter looking at collective modes and delta propagation in matter 
is presented in ref. J/J. A covariant formalism at finite temperature unifying 
the good features of the two formalisms is presented in ref. ||, with some 
application to nuclear matter within the relativistic Walecka model |J. 

Detailed analyses of nucleon properties along the lines of ref. ID] are avail- 



able at finite, but small temperatures T < 10 MeV [10], [TT|. Neither of the two 
approaches mentioned above is followed, but instead the smallness of the tem- 
perature justifies some approximations by means of which one finally follows 
the steps at T = substituting the Pauli distribution n(k) (1 for \k\ < k F , 
for \k\ > kp) by the thermal distribution 

n(k°) = [1 + expi^—^T 1 (1) 
k B T 

with ks the Boltzmann constant, k° = e(k) the nucleon energy and \x the 
chemical potential. The results of [|Hj are improved in by considering the 
correlation diagrams which lead to a large contribution to 7mE at momenta 
below the Fermi momentum. Similar steps are followed in ref. [[12] to deal 
with the propagation of the A at finite temperatures. 

The appeal of Thermo Field Dynamics is that one can continue to use 
Wick's theorem and the Feynman diagrammatic approach as at T = 0. The 
price is the introduction of an auxiliary field by means of which one constructs 
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a field doublet; the propagators become 2x2 matrices and the Feynman rules 
are now algebraic operational rules in the space of 2 x 2 matrices. 

In the present work we shall evaluate the imaginary part of the nucleon 
self-energy at finite temperature. In doing so we shall be following closely the 
steps of ref. where JmS was evaluated for nuclear matter at T = 0. 



The approach used in [O] was semiphenomenological, much in line with usual 



approximations done in the treatment of heavy ion collisions where the results 
obtained here are bound to be of relevance. The approach of [I3| evaluated 



correctly the second order nucleon self-energy diagrams but bypassed the use 
of an explicit NN potential. Instead, it used the fact that the sum of ladder 
diagrams in the low density limit provides the NN scattering t matrix, and 
its modulus squared, which appears in the evaluation of JmE, is related to 
the NN cross section. Hence the experimental cross section was used as input 
and the low density theorem [14], [15|] was automatically fulfilled. This theorem 
states that, as the nuclear density goes to zero, one has 

imEp^o = -„ a tot v p (2) 

where v is the velocity of the nucleon in the rest frame of the Fermi sea and 
(Jtot the total NN cross section. Long range correlations were also considered 
by means of an RPA sum with a phenomenological particle-hole interaction 
acting in the spin-isospin transverse channel. The approach proved to be nu- 
merically quite successful by comparing the results with those of elaborate 



and time consuming many body approaches like those of refs. [[U| PH . Spec- 
tral functions and occupation numbers were also evaluated in [jr3|] and were 
very similar to those obtained in other successful many body approaches with 
special emphasis in selfconsistency ]T$[. 

As with respect to the finite temperature treatment, we use the Thermo 
Field Dynamics approach by following the formalism, normalization and rules 



of ref. [19|, where a comprehensive and practical extract of this method is 
presented. 



2 Finite temperature formalism for the nu- 
cleon self-energy 

By following ref. 1TD| we have the thermal doublet for the nucleon field 

^ ( %) - I fH ) (3) 



i 



where ip(x) is the ordinary nucleon field and ip a support field for the formalism 
(t means the transposition with respect to the spinor index and f stands for 
conjugate). The thermal propagator is now defined as 

iG {a > b \x l7 x 2 ) =< 0(f3)\T[^ a \x l )^ b \x 2 )}\0([3) > (4) 
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where |0(/3) > is the state transformed from the vacuum by a Bogoliubov uni- 
tary transformation, and which has the property that the expectation value of 
an operator in this state is equivalent to its thermal average with the distri- 
bution of eq. (1). Hence the component of the Green's function, G^ u \ is the 
proper thermal average of the ordinary Green's function. 

The thermal free propagator, to be used in perturbation theory in the 
nonrelativistic approximation which we shall follow here, is given by 



with 



G^ ab \p) = G F (ab \p) + G° D {ab \p) 



(ab) , 



(5) 



qO (ab) 



(p) 



p° — e (p) + ^ e 



(6) 



GT\p) = ^(p°-e(p)) 



sin 2 9 p o ^sin20 p o 
^sin29 p o —sin 2 9 p o 



(7) 



cosd P° = (i + e -*)i/2 ; Sm9p0 



-x/2 



(l + e -x)i/2 



;x 



p — fJ> 

knT 



sin 2 9^ = = nip ) 



(9) 



This Green's function is in addition diagonal in spin but we omit the spin 
indices for simplicity. 

An alternative way of writing this propagator which is also suited to write 
the exact propagator is 







(ab) 



(10) 



with 



U F (p°) = 

and the exact propagator can be cast as 



cosOpO sinOpO 
—sinOpO cosOpQ 



(11) 



G ip) - U F (P ) ^ [p0 _ £{f) _ g^-i jU F ( P ) 

(12) 

with e(p) the kinetic energy of the particle and 

ReZ(p) = Re^(p) 

ImZ(p) = 7m£( n )(p)/cos20 p o 1 ' 
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In particular the G°^(p) component has an intuitive form given by 

r o(ii) / n = 1 - sin%o sin% 

yP) p° - e(p) - E p° - e(p) - E* 1 ; 

which in the limit of T = reproduces the standard form of the nucleon 
propagator in a Fermi sea. 

Hence in order to obtain the self-energy E(p) which renormalizes the nu- 
cleon propagator, only the thermal component E^ 11 ) needs to be evaluated. In 
the next section we show the approximation scheme that we follow to evaluate 
this magnitude. 



3 Semiphenomenological model for E 



This section follows closely ref. ||13|| . The generic Feynman diagram which 
we evaluate is the one in fig. 1, where the nucleon propagator in each of the 
baryonic lines is given by eq. (5). Note that in the limit of T = 0, and with 
the conventional separation of particles and holes, the usual polarization (fig. 
2a) and correlation (fig. 2b) graphs which lead to JmE JllJ are automatically 



generated (together with other two graphs with the interaction lines crossed 
which do not contribute to ME). 

In fig. 1 the indices a,b,c,d in the vertices are thermal indices. We are 
interested in E^ 11 ) and hence a = b = 1. Assuming for the moment the 
interaction lines to correspond to meson exchange, and considering also that 
no such mesons are present in the ground state of our many body fermionic 
system, the meson propagator would be diagonal in the thermal indices and 
hence c — d—1. 

We thus must evaluate the polarization function iT 11 ) (q) 

IL^(q) = -Aij ^ G° G°W(p + q) (15) 

where the factor 4 takes into account the sum over spin and isospin. 

Once again in the limit of T = 0, this polarization would account for the 
two diagrams in fig. 3, which are those accounted for by means of the Lindhard 
function j4j. 

By using G°( a ' b ) from eqs. (6,7) the p° integration can be easily performed 
and one obtains 



d 3 p j sin 2 d £{ jp ) cos 2 d e{ j; +ctj sin 2 9 e ^ +(f) cos 2 9, 



II (11) (g) = 4 / ^ <^ — IfMllf U£ (P +( i) _|_ 



e(p) 



(27r) 3 [q° + e(p) - e(p + q) + ie —q° - e(p) + e(p + q) + ie J 

(16) 

Next we evaluate E^ 11 ) corresponding to fig. 1 

Z^(k) =i\^ V\q) nM(?) G^\k - q) (17) 
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where V(q) would take into account the interaction due to the hypothetical 
meson exchange. 



Here again we follow the steps of ref. [13[] and sum the ladder diagrams 
which would replace V(q) by the scattering t matrix. (Note that medium cor- 
rections to t which would appear in the medium G-matrix are taken explicitly 
into account to second order with the structure of the diagram). We shall 



continue to use the same t matrix here. The studies of refs. [10], [TTJ, ^(J show 
indeed little dependence of the effective interaction on the temperature. 

In order to evaluate JmE 11 ' from eq. (17) a Wick rotation was made in ref. 
13fl which allows one to express /mE' 11 ' in terms of Imll^ u \ This is however 



not possible here because iT 11 ) from eq. (16) has overlapping cuts in the upper 
and lower half planes of the complex plane (unlike at T = where the cuts are 
confined to the second and fourth quadrant). However, an explicit evaluation 
of /mE' 11 ' is possible by first performing the q° integral in eq. (17) and then 
evaluating the imaginary part, with the result 

ImZW(k) = -47T / / (0| W° + sip) - e(k - q) - e(p + q)) 
■ { cos2e e(k^q) cos2 °e(p+q)Sin 2 e £{ ^ - sin 2 6 £{l :_ (t) sin 2 e £ ^ +t f ) cos 2 e £{ ^} 

The spin-isospin averaged value of \t\ 2 assumed in eq. (18) is then replaced 
by 4ti<tnn/M 2 , where M is the nucleon mass and <jnn the spin-isospin averaged 
NN elastic cross section. Since pion production is not explicitly taken into 
account, this restricts us below the pion production threshold. The final step 



in ref. Jl3fl is to consider the polarization or RPA corrections to the interaction. 

The consideration of the polarization was an important ingredient in ref. 
13fl , which reduced /mE particularly at small energies, and provided results 



similar to those found in refs. |16|, [L7|]. We implement it here too. For this 
purpose we need to evaluate H^ u '(q), both the real and imaginary part, which 
cares about ph excitation, and n^ 1 ^), the equivalent term accounting for 
Ah excitation. At T = these quantities are the familiar Lindhard functions 



Un(q), U^iq), respectively, used in ref. |L3 



The real part of Un{c[), unlike ImllN(q), is not affected by Pauli blocking 
0], hence finite temperature, which modifies occupation numbers, has not 
much of a consequence in the change of ReU]\r(q). On the other hand there is 
no Pauli blocking in the Ah excitation since we do not have a Fermi sea of A's. 
For these reasons we keep iten^ 11 ) and at finite temperatures equal to 
ReUN^.U^q) at zero temperature. However, we evaluate Mn' 11 ^) from 
eq. (16) at finite temperature. The reason is that keeping /mil' 11 ' ^ is 
important in order to avoid singularities coming from poles of zero sound (q° 
proportional to |g| at small energies) which are strongly dumped at finite T. 

The expression for ImUS 11 ^ obtained from eq. (16) is given by 

ImU^\q\ q) = -\ dp^{sin 2 6 e{f) + sin 2 e em+q » 

(19) 

-2sm 2 6' e(??) sm 2 6' e ( ??)+? o} 



6 



where 
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Pmin= -\q-^-\ (20) 

p 2m 

The polarization correction replaces the interaction by the induced interac- 
tion [21] (see eq. (15) of ref. |13|]). Furthermore, we can perform some trivial 



integrals and eliminate the 5 function with the final result 
ImEV\k) = -fp / °° Qdqf-i dcos6 f °° pdpQ(l - A 2 ) 

' \l-Vt(q)U(q)\* \q°=k°-e(k-q) (21) 

' { cos2 ^e{k~q) cos2 ^e(p+q)Sin 2 e £ ^ - sin 2 6 £ ^_^sin 2 6 £ ^ +(i) cos 2 6 £ ^} 

where Vt(q) is the transverse part of the spin-isospin interaction and U(q) = 
n^ 11 )(g) + Il^^g). The arguments leading to this modifications and expres- 



sions for V t (q) and £/jv(g), U&{q) can be found in ref. and we do not repeat 
them here. Furthermore the angle in the integral over cosO in eq. (21) is the 
angle between q and k. The magnitude A in the argument of the step function 
is given by 

f + kqcos9) (22) 




2M M M J 

with k, q the modulus of k and q respectively. 

The value of the chemical potential function of the density and T is 

obtained, as usually done, by the normalization condition 

r d 3 k 1 
P = A J (2vr) 3 1 + exp[(e(k) - fi)/k B T] ^ 

Eq. (21) provides TmE^ 11 ) as a function of k°, k. In the next section we 
show the results which we obtain for this quantity. 



4 Results and discussion 

In fig. 4 we show the results of — JmE at T = for two densities, p 
(0.17/m~ 3 ) and po/2, obtained with the present formalism in the limit of 



T = 0. The results agree with those in ref. |13] calculated with the T = 
formalism and also with those of the microscopic approach of ref. |IE |. Note 
that since only kinetic energies are used as input to evaluate X^ 11 ), the value 
of n is referred to an origin of energies at \k\ = 0. We are not interested in 
i?eS and to overcome the arbitrary origin of energies we plot the magnitudes 
in terms of u — fi ( u = k° ). In fig. 4, \k\ is taken at the value y2Muj. 



This justifies small differences with ref. [13| where the value of \k\ satisfying 



the dispersion relation between \k\ and u was taken. In fig. 4 we observe the 
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typical features that JmS is proportional to (u — p) 2 . In the calculations we 
find that /mS' 11 ' changes sign at u — p,, with /mS' 11 ' < for uj > p. In this 

case 




(24) 



cos29 k o = 1 — 2sin 2 6 k o 

= 1 - 2n(k°) T=0 = 

and then 

imE 11 

Im£ = — — = -\ImE n \ (25) 

cos 2 9 k o 

The Green's function G^ u \ by using eq. (12), will then be 



Q(uj — p) 6(/i — uj) 



k° - e{k) - i?eS( n )(A;) + i|/mE( n )| k° - e{k) - ReE( n ){k) - i|im£ (11) | 

= i (26) 

k° - e{k) - S( n ) 

as it should be. 

In fig. 5 we plot JmS' 11 ' at p = po as a function of uo — fi, with |fc| = \j2Muj 
and /i calculated from eq. (23). As can be seen in the figure, JmS^ 11 -' is 
always zero at u — [i. However, ImS(//, k) is different from zero at finite 
temperatures, contrary to the situation at T = where it is zero. In order 
to envisage this we see that the evaluation of JmS from eq. (13) involves a 
fraction of the type 0/0 which we determine using l'Hopital rule and find 



/m2(„, *) = lim ^"^ff = k)\ k ^ (27) 

We can see that JmS(/i, k) vanishes at T = 0, as we already said. 

In fig. 6,7,8 we show the results for -ImE as a function of u — fi for 
different temperatures, calculated for three different densities, po/2,p and 
2po- We can appreciate that as T increases -/mil also increases in all the 
range of energies calculated there. We should note that in evaluating JmE, 
the factor cos29 k a appearing in the denominator of eq. (13) is very important 
and makes imS ^ at T ^ for all the range of energies, while JmE^ 11 ' passed 
through zero. This is a genuine temperature dependent property which would 
be lost if a T = formalism, changing the Fermi distribution by the thermal 
distribution of eq. (1), were used. 

At this point it is interesting to compare our results for JmE with those 
which we would obtain using standard formulae of Fermi-liquid theory 



The formula used there in our notation for JmE, removing the cut off in the 
integral, would be given by eq. (18) changing the minus sign in the curled 
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bracket {cos 2 and sin 2 terms) by a positive sign. Instead our formula for 
ImE uses eq. (18) with a minus sign (which provides ImE^) and then we 
divide by cos26 po (as shown in eq. (13)) in order to obtain /mS. The same 
prescription would be taken if one uses instead eq. (21) which incorporates 
the effects of polarization. 

It is easy to see that at T = both formulae are identical. Indeed for 
u > fi only the first term in the curled bracket of eq. (18) (or (21)) (the one 
with cos 2 cos 2 sin 2 ) contributes, while for u < \x only the second term in the 
bracket (the one with sin 2 sin 2 cos 2 ) contributes. Then when dividing /mS' 11 ' 
by cos29 p o, given by eq. (24), we obtain a formula for ImE given by eq. (18) 
where the minus sign in the curled bracket is changed to a positive sign, exactly 
the formula used in Fermi-liquid theory [p2]| . 

The identity of the two formulae holds, however, at any temperature. This 
can be seen analytically using eqs. (8) for sin8, cos9 and x, hence 



(cos 2 9 £l cos 2 9 £2 sin 2 9 £3 — sin 2 9 £l sin 2 9 £2 cos 2 9 £3 ) N 

o\k + e 3 — e\ — e 2/ 



cos29 k o 

1 e ^3^ e (xi+X2-X 3 ) _ e X k0 _|_ J 

k B T (1 + e Xl )(l + e X2 )(l + e x *) e x ^ - 1 
1 e Xl e X2 + e x:i 



5(x k o +x 3 —xi — x 2 ) 



7 7 /1 .Ml \fl -~^ S ( X kO+X 3 -X 1 -X 2 ) 

KbT (1 + e Xl )(l + e X2 )(l + e X3 ) 
= (cos 2 9 £l cos 2 9 £r sin 2 9 £3 + sin 2 9 £l sin 2 9 £2 cos 2 9 £3 )5(k° + £3 — e\ — e 2 ) 

(28) 

where the constraints of the ^-function have been used in the second step. 
This is an interesting finding which stresses the value of the results obtained 
in Fermi-liquid theory based on the concept of quasiparticles, by comparison 
to a method in principle more microscopic, like the one used here. 

Next we would like to extract a practical magnitude from these results 
which can be used in calculations of heavy ion collisions. Recall that in the 
semiclassical approach one has 

ImE = — — (Jnn v p (29) 

One can then define an effective NN nucleon-nucleon cross section at different 
T and p by means of 

pf f ImE , , 

af N = -2 (30) 

vp 



as done in ref. [^3|], where p defines a probability of collision per unit 
length for the nucleon. In order to facilitate the use of this magnitude we 
have parameterized va^N, with v = \k\/M for the different densities and 
temperatures evaluated here. We take the following functional form 
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vo e M = Y,a n {p,T)u n (31) 

n=0 

where uo is the nucleon kinetic energy, k 2 /2M . The fit is valid for values of 
uj > p in figures 6,7,8. The coefficients a n (p, T) are given in tables I, II, III. One 
can obtain va^N f° r intermediate values of p and T interpolating between the 
values of va e ^ given by eq. (29). 

5 Conclusions 

We have used a model to evaluate JroS for a nucleon in nuclear matter at finite 
temperatures, which combines the formalism of Thermo Field Dynamics with 
empirical magnitudes of the NN interaction. This model at T = coincides 
with a semiphenomenological approach studied earlier, which proved rather 
successful in reproducing nucleon properties in nuclear matter obtained with 
more microscopical approaches. 

We have obtained ImT, for different values of the nuclear density and sev- 
eral temperatures. One of the interesting findings is that JmS grows steadily 
with the temperature. The changes produced by the temperature are more 
striking at energies around the chemical potential where ImE is zero at T = 
and takes finite values at T ^ 0. 

We found that the genuine effects of the temperature, given naturally in 
the formalism of Thermo Field Dynamics, were important, and the differences 
with respect to simple calculations, where n{k) at T = is replaced by the 
thermal distribution, can be appreciable. In order to facilitate the use of the 
results obtained here, we have parameterized them by means of easy analytical 
formulae. The parameterization is given for an effective NN cross section in 
the medium, such that p gives the probability of collision per unit length 
for a nucleon in the nuclear medium. With the given formulae one can easily 
interpolate the results and obtain this magnitude for different densities and 
temperatures. These results should be useful in the analysis of heavy ion 
reactions. 
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6 Table Captions 



Table I Parameters of eq. (28) to evaluate va e ^ } N at p = po/2. The parameters 
a n have dimensions such that, with uo given in MeV, the results for va e J^ are 
given in mb. 

Table II Same as table I for p = p . 
Table III Same as table I for p = 2p . 
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7 Figure Captions 



Fig. 1. Generic Feynman diagram to evaluate the nucleon self-energy. The 
indices a, b, c, d are thermal indices. The nucleon propagator corresponding to 
the baryonic lines is given in eqs. (5-9). 

Fig. 2. a) polarization, b) correlation graphs contributing to ImT, at T = 
and contained in fig. 1. Here the direction of the arrows stands for the 
conventional hole ( down ) and particle ( up ) propagators. 

Fig. 3. Polarization graphs appearing at T = with the same notation for 
the lines as in fig. 2. 

Fig. 4. —ImE(ui,k) at T = as a function of uj — p, with k = V2Mu, 
evaluated for two densities. 

Fig. 5. ImE' n '(w, fc) at p = p for several temperatures as a function of 
uj — p with k = V2Mu. The solid line is for T = 0. The other curves 
correspond to T = 2 MeV ( long dashed-dotted line ), T — AMeV ( dashed 
line ), T = lOMeV ( dotted line ) and T = 20 MeV ( short dashed-dotted 
line ). At values uj — p < 0, they appear correlatively with increasing values 
of JmE' 11 ' as T increases. 

Fig. 6. -/mE(w, k) at p — po/2 as a function of uj — p for k = \j2Muj for 
several temperatures T = 0, 2, 4, 10, 20 MeV with the same notation as in fig. 
5. —ImT, increases with increasing T. 

Fig. 7. Same as fig. 6 at p = Pq. 

Fig. 8. Same as fig. 6 at p = 2p . 
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Table I. 



T 

a n 


MeV 


2 MeV 


4 MeV 


10 MeV 


20 MeV 


a Q 


-2.549 


-2.305 


-1.425 


0.982 


5.107 




9.658 10- 2 


9.098 10- 2 


7.043 10- 2 


3.404 10~ 2 


-2.170 10~ 2 


a 2 


-9.706 10- 5 


-3.996 10- 5 


1.569 10~ 4 


4.367 10~ 4 


8.413 10~ 4 


a 3 


-4.349 10~ 7 


-6.674 10~ 7 


-1.435 10~ 6 


-2.408 10- 6 


-3.857 10- 6 


a± 


1.135 10" 9 


1.456 10~ 9 


2.487 10^ 9 


3.724 10^ 9 


5.709 10~ 9 



Table II. 



T 

Qj n 


MeV 


2 MeV 


4 MeV 


10 MeV 


20 MeV 


a 


0.553 


0.411 


0.793 


1.695 


3.192 


ai 


-3.541 10~ 2 


-2.722 10~ 2 


-3.709 10~ 2 


-4.496 10~ 2 


-4.430 10~ 2 


a 2 


7.790 10" 4 


6.834 10~ 4 


7.846 10" 4 


8.343 10~ 4 


7.542 10" 4 


a 3 


-2.860 10~ 6 


-2.458 10~ 6 


-2.861 10~ 6 


-3.027 lO" 6 


-2.651 10~ 6 


04 


3.521 10~ 9 


2.963 10~ 9 


3.501 10~ 9 


3.719 10~ 9 


3.215 10~ 9 



Table III. 



T 

(l n 


MeV 


2 MeV 


4 MeV 


10 MeV 


20 MeV 


a 


1.493 


1.382 


1.347 


2.205 


2.861 


a\ 


-5.370 10~ 2 


-5.017 10~ 2 


-4.795 10~ 2 


-6.283 10~ 2 


-6.065 10~ 2 


a 2 


5.866 10~ 4 


5.555 10~ 4 


5.355 10~ 4 


6.533 10~ 4 


6.214 10" 4 


a 3 


-1.659 10~ 6 


-1.553 10~ 6 


-1.484 10~ 6 


-1.874 10~ 6 


-1.770 10- 6 


CI4 


1.695 10- 9 


1.570 10~ 9 


1.488 10- 9 


1.940 10~ 9 


1.834 10~ 9 
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Fig. 6 




(o-|i (MeV) 

Fig. 8 



Fig 1 




